Field-Induced Degeneracy Regimes in Quantum Plasmas 

M. Akbari-Moghanjoughi 

Azarbaijan University of Tarbiat Moallem, Faculty of Sciences, 
Department of Physics, 51745-406, Tabriz, Iran 
(Dated: August 28, 2012) 

Abstract 

It is shown that in degenerate magnetized Fermi-Dirac plasma where the electron-orbital are 
quantized distinct quantum hydrodynamic (QHD) limits exist in which the nonlinear density waves 
behave differently. The Coulomb interaction among degenerate electrons affect the electrostatic 
nonlinear wave dynamics more significant in the ground-state Landau quantization or the so- 
called quantum-limit (/ = 0) rather than in the classical-limit (/ = oo). It is also remarked that 
the effective electron quantum potential unlike the number-density and degeneracy pressure is 
independent of the applied magnetic field in the classical-limit plasma, while, it depends strongly 
on the field strength in the quantum-limit. Current findings are equally important in the study of 
wave dynamics in arbitrarily-high magnetized astrophysical and laboratory dense plasmas. 

PACS numbers: 52.30.Ex, 52.35.-g, 52.35.Fp, 52.35.Mw 
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I. THEORETICAL REVIEW 



Plasmas under extreme conditions such as high pressure and magnetic field has well- 
known applications in astrophysics [l-4|. Since the nobel work of Chandrasekhar, it is well 
understood that the thermodynamic properties of degenerate matter can be fundamentally 



altered due to relativistic consideration caused by gigantic gravity force 



a. 



.t has been 



shown that the relativistic degeneracy can affect the whole thermodynamics [6| as well as 
linear and nonlinear wave dynamics [7H9| of Fermi-Dirac plasmas. In fact the remarkable 
property of relativistic degeneracy has been shown to be the lead cause of complete chain 
of stellar evolution and existence of variety of compact stars in our Universe 10|]. The 



associated magnetic field with some highly dense stellar objects 
star forming molecular clouds 112 



11| its inevitable role in 



- |l4| has brought the attention to the problems such as 



behavior of atomic structure and Fermi-electrons under an arbitrary strength magnetic field. 



The earliest extensive study of this kind belongs to Canuto and Chiu 
shown 



17 



15 



|22|. It has been 



23l | that the electronic density of state, hence all thermodynamics quantities, 



become quantized due to the familiar Landau orbital quantization. 

Obviously, the significant change in equation of state of ordinary matter has yet to be 
observed in laboratory due to high magnetic fields close to critical value. Be = mlc^/eh ~ 



24j of quantization given rise to De 



4.41 X lO^^G, required, however, the oscillatory nature 
Haasvan Alphen effect in many experiments is an undeniable proof of the Landau theory, a 
complete review of the properties of degenerate matter under arbitrary magnetic field may be 



found in literature l|, l25[. One of the remarkable predictions based on the Landau extension 



of Fermi-Dirac statistics is the existence a metastable ferromagnetic state called LOFER 



26- 



30| which might explain the unsolved origin of high magnetic field in dense astrophysical 



plasmas 27|, |30(|. The studies based on LOFER condition estimates magnetic field of order 



10**G for white dwarf and of lO^'^G for neutron star electron-number density ranges 31 |. 
Based on the magnetohydrodynamic formalism it has been recently shown that the LOFER 
concept for the astrophysical density ranges can lead to quantum collapse where the Fermi- 
Dirac degeneracy pressure is totaly canceled by the Landau spin-orbit magnetization pressure 



32 



33| . Future progress in intense laser technology may facilitate the observation of Landau 



quantization in strong laser-matter interaction experiments 



351 ]. However it is thankful 



of nature for providing us with the astrophysical laboratory to observe phenomena in which 
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quantum, relativity and magnetism effects are all involved. 

In current research we examine the effect of such extreme conditions on nonlinear wave 
dynamics incorporating the relativistic Coulomb interactions which is most relevant to highly 
degenerate plasmas such as the crystalized core of white dwarfs or neutron stars. It has been 
shown that the magnetic field and Coulomb interactions may lead to minor deviations from 
the Chandrasekhar mass-radius relation [sgI. There are some other effects such as electron 
distribution non-uniformity, exchange and correlation effects which are of less importance 
compared to the Coulomb interaction [37] and are neglected in the present work. The pre- 
sentation of the article is as follows. The basic normalized plasma equations are introduced 
in Sec. HIl A general nonlinear arbitrary-amplitude solitary solution is given in Sec. IIIII 
The numerical presentation and discussion is given in Sec. |IV]and conclusions are drawn in 
section |Vl 



II. QUANTIZED FERMI-DIRAC PLASMA MODEL 

In this section we present simple quasineutral quantum hydrodynamics model for a com- 
pletely degenerate quantized Fermi-Dirac plasma of with interacting electrons. For simplic- 
ity, we consider only perturbations parallel to an arbitrary-strength uniform magnetic field 
applied in the z direction to probe the possible magnetic quantization effect on nonlinear 
density structures propagating in z direction parallel to the existing field. Based on the 
standard quantum hydrodynamics (QHD) model, the continuity and momentum equations 
for electrons and ions ignoring the Bohm force on ions, can be combined to give a complete 
set of equations in the plasma center of mass frame reference, as 

| + V-(pu) = 0, (1) 

where, p = {mini + rriene) ~ mjnj and, u = {nerrigUe + nirriiUi)/ p are the center of mass 
density and velocity of quantum plasma, and 

dn ^ dPtot ^ pff d r 1 d^^A 
^ dt dx 2memi dx \ ^/n^ dx"^ J ' 

where, Rtot = Pe\\ + Pi + Pint with P^w, Pi, Pint being the parallel component of the field- 
dependent degeneracy pressure of electrons, classical ion pressure and interaction pressure 
such as Coulomb, respectively and h is the scaled Plank constant. It is known that in the 
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Fermi-Dirac plasma among the wide variety of interactions the Coulomb interaction is the 
most dominant [36] which will be considered in our problem. We further use the following 
scaling to reduce to normalized model 



X — )• X, t — )■ , p — )■ PPq, U — )■ UCs, 



where the normalizing factors, po, u. 



(3) 



pj — e/mi^JA^lpQ and Cg = c-^me/rrii denote the equi- 
librium plasma mass-density, characteristic ion plasma-frequency and ion quantum sound- 
speed, respectively. The bar notations denote the dimensionless quantities, hence, ignored in 
forthcoming algebra for simplicity. The normalization of the continuity equation is straight. 
On the other hand, the quasineutrality condition (ng ~ = n) along with the scalings 
introduced in Eq. ([3]) gives rise to the following equation for the momentum 

dt 



(4) 



p dx ' dx \^/p dx"^ 
where H = ^Jmi/2m^{hJJpi) / {rrieC^) is the relativistic quantum diffraction parameter. The 
thermodynamic quantities of a Fermi-Dirac gas in an arbitrary strength magnetic-field is 
given by Canuto and Chiu 17|, |36|]. The quantized fermion number-density (not normalized) 



reads as 



neiepcl) = nc7_^ (2 - 5ifi)^J e],^ - 2/7 - 1, (5) 
z=o 

where, ric = mlc^ / 27r'^ Jr' , = Bq/ Be with Be = mlc^/eh ~ 4.41 x lO^^G being the fractional 



critical-field parameter, epe = Vl + -R^ = Epe/^eC^ is the normalized Fermi-energy and 
R = PFe/^eC is the normalized Fermi-momentum the so-called relativity parameter. Also, 
Im is the maximum filled Landau-level defined as Im = {^'pe ~ l)/27 > /• The parameter 61^ 
is one for the quantum-limit (/ = 0) and zero for other Z-values. On the other hand, the (not 
normalized) quantized pressure can be written as 



Pe\\ {SPe, 7) = |nc7meC2 E (2 - 5i,o) \ SFe^/s^e "2/7-1 

1=0 

£Fe+A/£L-2^7-l 



;i + 2/7) In 



(6) 



The small-field expansion for the critical-field parameter, 7, can be written as 

rie(i?,7) = K + g + 0(7^)) , 
Pe\i{R,-f) = I + i?2(2i?2 - 3) + Bsinh-^i?] 

+ 0(7^ 



(7) 



R 



2sinh^^i? - 1 + 



It is noted that the 7 = hmit reduces to the Chandrasekhar equation of state. 



23|. 



The above equation of state may also be written in terms of Hurwitz zeta-functions 
For instance, for the particle density and the pressure of Landau quantized Fermi-Dirac gas, 
we may write 



n,(i?,7)=n,(27)3/2i/-i/2 

Pell (i?, 7) = ^(27)^/^0^ ^rf?, .g. 

H,{q) = h{z,{q})-h{z,q + l)-lq~\ 

00 

h{z,(l) = E {n + qy. 

n=0 

where h{z, {q}) is the Hurwitz zeta-function of order z with the fractional part of q as 



argument. It has been shown [23[ that the thermodynamics parameters may be separated 
into monotonic and oscillatory parts. Figure 1 (left plot) shows the exact quantized particle- 
density (thick curve) along with the monotonic and oscillatory components for 7 = 0.5. Also, 
Fig. 1 (right plot) reveals that up-to i? = 1 (/ = 0) the dependence is linear and the step-like 
feature introduced in density structure due to the Landau quantization is smeared-out in 
the classical- limit, / = 00. The expansion of quantized density given in Eq. ([8]) for the 
classical- limit, -R^/27 ^ 1, |23| leads to 



2 

neiR,l) = -Uc 



+ 7V4i? + l{2jf'h{-l/2, {PV27}) + 0(7^) 



(9) 



The first two terms correspond to the monotonic (shown in Eq. ([7])) and the third term 
belongs to the oscillatory one. There are distinct critical-field regimes which are apparent in 
Fig. 1, namely, the quantum-limit, / = (i.e. -R^/27 ^ 1) and the classical-limit, / = 00 (i.e. 
R^/2'~f ^ 1). These regimes for 7 = 0.5 case are denoted as i? < 1 for the quantum-limit 
and -R > 3 for the Chandrasekhar classical-limit, in Fig. 1. The magnetic-field and density 
range shown in both limits in Fig. 1 reveal the importance of both in the astrophysical and 
laboratory applications. By defining a dimensionless normalizing parameter Rn = {no/ncY^^ 



we have for the classical Chandrasekhar-limit (/ = 00), n'^R^ = R^ 15|, |3J] and for the 
quantum-limit, we obtain, n'^Rl = 7-R (e.g. see Eq. ([5]) for / = 0). 

We now return to the normalized hydrodynamic equation in terms of effective potential 
instead of pressure. Ignoring the ion pressure compared to that of degenerate electrons, we 
have 

^ = + 4„ = + (10) 

at ox ox VyP / 



where, the normahzed relativity parameter Rq and the relativistic diffraction param- 
eter, H = \/mi/ 2me{hoOpi) / (meC^) , are related through the simple relation H = 



Then in normalized form we have 



rae(i?,7) 



d_ 



H-i/2{q) 



dq 



dr] 
dR 



dR, r] 



27' 



(11) 



This gives a remarkably simple form as $e = ^/l + R^ = Spe for any quantization level. This 



remarkable feature indicates that in the classical-limit {n = nf = [R/Rq)^) the effective 



potential is completely independent of the magnetic field strength and only depends on the 



fermion number-density. However, in the quantum-limit {n 



potential depends strictly on the fractional field-parameter, 7, i.e. 



where, we have redefined Rq = {no/ncY^^ = (po/ PcY^^ with pc = m. 




n 



On the other hand, for the interacting Coulomb pressure we have 36 1 



jR/Rl) the effective 
(12) 

; ~ 1.98 X lO^gr/cm^. 



Pc{R.i) 



,5^2/3 



10L4 



a 



37r 



nl/3 



(13) 



rae(i?,7)_ 

where, the parameters, a = e^/hc ~ 1/137 and Z are the fine-structure constant and the 
atomic number, respectively. Hence, we can write for both / = 00 and / = regimes 

1 dPc{R,^) 



int 



-dR = -/3Rop 



1/3 



(14) 



n(i?,7) OR 

where, /3 = a{2^^^/5){3Z'^/7iY^^. Note that this effective potential introduces a negative 
force in the momentum equation unlike the electron degeneracy pressure. 



III. SOLITARY EXCITATIONS IN QUANTIZED DEGENERACY REGIMES 

In order to obtain a stationary soliton-like solution we follow the conventional pseudopo- 
tential approach. In the moving frame with speed M we take the coordinate transformation 
C, = X — Mt which brings us to the co-moving soliton frame. Introducing the new coordi- 
nates into the normalized hydrodynamic equations, integrating by use of the appropriate 
boundary conditions, lim p = 1 and lim = 0, we solve the continuity relation to give 
u = M [l/p — 1). Eventually, by making use of this relation in the momentum equation and 
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defining the new variable, p = A^, we arrive at the following single differential equation for 
the classical- limit, we have 



(15) 



and, for quantum-limit we may write 

= — (1 - -M\l- A-') 



1 + 



l + ^-f3RoA^/^+/3I^. (16) 



^ V T 

Now multiplying of the Eq. f[T^ and Eq. f lTC]) by the quantity dA/dS, and integrating with 
respect to, ^ with the mentioned boundary conditions, lead us to the well-known energy 
integral of the forms given below in terms of plasma mass-density 



{d^pf/2 + U{p) = 0. 
Hence, for the classical- limit, we have 

UM = {^M^RliP - 1)' + R^P [3 [V^TM - sji + i?gp2/3pi/3^ 

+2RI [v/rT^(4p - 1) - 3pv/l + RIp^'"" + /3^o(l - 4p + Sp^/^)]" 
-3psinh'^i?o + 3psinh"^(i?oP^/^)} , 
and, for the quantum-limit it follows that 

2M2(1 - pf - 27- VVt' + RIp" + 4pVV7' + ^oP' 



(17) 



(18) 



f/o(p) 



1 I p3 

2^2 1 -^0 



+i?o/3p(l - 4p + 3p4/3) - 2pVV7' + ^oP' 
-F27P [sinh"^(7~ii?ij) - sinh'^ {p^i'^ Rl)] ] . 



(19) 



Existence of the localized density profiles requires the following basic conditions to met all 
together 



U{p)\ 



dU{p) 



p=i 



dp 



0, 



d^U{p) 



p=i 



dp^ 



< 0. 



(20) 



p=i 



Evidently, the first two requirements are satisfied for the pseudopotential given by Eqs. fllSp 
and f|T9|) . The third condition, however, can be directly evaluated using the pseudopotentials 
for both quantized degeneracy. The soliton, therefore, exists if the there is pseudopotential 
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root other than at the unstable point p = 1 which will be examined later. In such conditions, 
the spatial extension of the solitary wave in then given as 
dp 



± 



(21) 



'1 V-2t^(p) 

we now return to the evaluation of the second derivative of the Sagdeev potential in both 
degeneracy limits, at unstable density p = 1, which leads to the following expression for the 
classical-limit 



p=i 



+ Ro\ (3 



Ro 



(22) 



and the following relation for quantum-limit 



d'^Uoip) 



dp^ 



3H^ 



+ Ro\ (3 



3Rl 



(23) 



The soliton Mach-range is therefore defined through the following critical Mach-values in 
each degeneracy limit. For classical-limit, we have 



M 



R2 



RoP 



s^iTM 3 • 

Note that, the given relation reduces appropriately to the one in Ref. 
(/3 = 0) limit. On the other hand, for the quantum- limit, we obtain 



(24) 



38 



39| in the field-free 



M' 



Rl 



RoP 



(25) 



i^^F+M 3 • 

Now let us examine the possibility of a root other than p = 1 which is other essential 
condition for existence of potential valley. For the classical-limit, it is found that 



M2 

limf/oo(p) = -77 > 0, hm f/oo(p) 



-00 X sgn(l — f3), 



(26) 



which indicates that a rarefactive solitary excitation always exists with Mach-values below 
that of critical value defined above, and, for the quantum-limit, we have 



M2 

lim[/o(p) = — > 0, lim Uoip) 

p— >-0 M p— s>oo 



-00, 



(27) 



revealing also the fact that only rarefactive density profile is possible in this limit. In the 
next section we bring into the attention the possible differences in soliton characteristics 
caused by distinct magnetic degeneracy limits and discuss the variations of solitary profiles 
in terms of change in various fractional plasma parameters. 
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IV. NUMERICAL INSPECTION 



In this section we examine the characteristics of sohton dynamics in the two magnetic 
degeneracy hmits under the influence of change in many plasma parameters such as fractional 
field-parameter, 7, the normalized relativity parameter, Rq, and the atomic- number, Z. 
Figure 2 shows the soliton stability volume in classical-limit, i.e. when a large number of 
Landau orbital-lcvcls arc filled. It is observed that the increase in the fractional Fermi- 
momentum Ro widens the soliton Mach-range, while, the change in the atomic number does 
not affect the soliton stability significantly, in this regime. Comparing this figure with Fig. 3 
which is the corresponding volume in the case of quantum-limit {I — 0), reveals that, in this 
case the atomic number has significant effect on the soliton stability. In fact it is observed 
that, for this regime at low relativity parameter region, there are unstable atomic-number 
values for a given value of Rq. It is remarked that for / = limit the less the atomic-number 
the more stable density structures are. 

Figs. 4 and 5 depict the variations in pseudopotentials for the cases of Z = 00 and I — 
while one of the plasma parameters is changed and the others are fixed. The comparison 
of the figures reveal the fundamental differences on plasma parameters on potential profiles 
in the two limits. Both potential width and depth are affected relatively strongly in the 
quantum-limit degeneracy. The corresponding rarcfactive solitons shown in Figs. 6 and 7 
with identical parameter values to those in Figs. 4 and 5 for two regimes confirm the above 
statement, clearly. While the effect of the plasma parameters on soliton profiles are similar 
in the two cases, the nature of the effects are fundamentally different. For I — 00 the increase 
in atomic-number (Fig. 6(a)), relativity parameter (Fig. 6(b)) and fractional field-strength 
(Fig. 6(c)) all widen the soliton width without altering the amplitude significantly. Also, 
increase in the Mach- number decreases/increases the soliton amplitude/ width slightly for 
this case (e.g. see Fig. 6(d)). 

On the other hand, for the case of i = 0, the increase in the atomic-number decreases 
the soliton amplitude (Fig. 7(a)), while the increase in the normalized Fermi- momentum 
increases the soliton amplitude (Fig. 7(b)). Also, it is observe from Fig. 7(c) that the in- 
crease in fractional field strength decrease the sohton amplitude significantly. Moreover, the 
effect of increase in soliton Mach value (Fig. 7(d)) strongly decreases the soliton amplitude 
in this case. All these characteristics differ fundamentally in the two magnetic quantization 
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limits. These findings are quite comparable to the distinctive features of nonlinear wave 
dynamics in nonrelativistic and relativistic degeneracy limits reported previously for the 
field- free Fermi-Dirac plasmas 38 1. 



V. CONCLUSIONS 



Based on quantum hydrodynamics formulation and using pseudopotential approach it was 
shown that for a Coulomb interacting Fermi-Dirac plasma under arbitrary strength uniform 
magnetic field two distinct degeneracy limits exist in which the parallel electrostatic solitary 
wave dynamics are fundamentally different. It was further revealed that plasma parameter 
such as atomic-number, fractional Fermi-momentum, fractional field-strength and soliton 
Mach affect the wave characteristics significantly only in the quantum-limit degeneracy 
compared to that of classical-limit. These findings can help better understand the physical 
processes in astrophysical compact objects as well as the strong laser-matter interactions. 
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FIGURE CAPTIONS 



Figure- 1 

The monotonic and oscillatory components in fermion number-density and the exact den- 
sity (thick curve) is shown in left plot, indicating a linear region of dependence on fractional 
electron relativistic Fermi-momentum, R, {ug oc jR) which corresponds to quantum-limit 
Landau-quantization, I — 0. The plot to the right shows the fading of the oscillatory part and 
asymptotic Chandrasekhar-limit approximation (ug oc R^) when the Fermi-Dirac plasma is 
in classical Landau quantization-limit (/ = oo). 

Figure-2 

Figure 2 shows a volume in M-Z-Rq space for classical-hmit, {I — oo), in which a localized 
magnetosonic density localized excitation can exist. 

Figure-3 

Figure 3 shows a volume in M-Z-Rq space for quantum-limit, {I 
magnetosonic density localized excitation can exist. 
Figure-4 

The variations of Sagdeev pseudopotentials of localized density-structure for the classical- 
limit magnetic Fermi-Dirac plasma with respect to change in each of several independent 
plasma fractional parameters, namely, normalized soliton-spccd, M, the atomic-number, Z, 
fractional field-strength, 7, and the normalized relativity parameter, Rq, while the other 
three parameters are fixed. The dash-size of the profiles increase according to increase in 
the varied parameter. 

Figure-5 

The variations of Sagdeev pseudopotentials of localized density-structure for the 
quantum-limit Fermi-Dirac plasma with respect to change in each of several independent 
plasma fractional parameters, namely, normalized soliton-speed, M, the atomic-number, Z, 
fractional field-strength, 7, and the normalized relativity parameter, Rq, while the other 



= 0) , in which a localized 
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three parameters are fixed. The dash-size of the profiles increase according to increase in 
the varied parameter. 

Figure- 6 

The variations of rarefactive localized density-structure shape for the the classical-limit 
magnetic Fermi-Dirac plasma with respect to change in each of several independent plasma 
fractional parameters, namely, normalized soliton-speed, M, the atomic-number, Z, frac- 
tional field-strength, 7, and the normalized relativity parameter, Rq, while the other three 

parameters are fixed. Identical values as Fig. 4 is used in this figure. The thickness of the 
profiles increase according to increase in the varied parameter. 

Figure- 7 

The variations of rarefactive localized density-structure shape for the quantum-limit mag- 
netic Fermi-Dirac plasma with respect to change in each of four independent plasma frac- 
tional parameters, namely, normalized soliton-speed, M, the atomic-number, Z, fractional 
field-strength, 7, and the normalized relativity parameter, Rq, while the other three param- 
eters are fixed. Identical values as Fig. 5 is used in this figure. The thickness of the profiles 
increase according to increase in the varied parameter. 
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